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verview  and  motivation 

UCAV  Simulation  Issues 
Simulation  hierarchies 


♦  Static  Case  Validation  of  DES 

♦  Forced  Motion  Validation  of  DES 


♦  Embedded  LES  Modifications  to  DES 


Future  Areas  of  Research  Necessary 


♦  Conclusions 


UCAV  Simulatiofffil 
Issues 

>  Unmanned  Combat  Air  Vehicles  are  capable  of  supeg-maneuverabilty 

>  Main  Challenges 

>Maneuvers  occur  at  high  Reynolds  numbers  for  which  the  underlying  fluid 
motion  is  usually  turbulent 

> Incorporates  massively  separated  flows  and  complicated  vortical  flows 
>Complete  simulation  requires  solid-body  motion|6-DOF,  and  aeroelasticity 

>  Wind  tunnel  tests  problematic 

>  Important  Reynolds  number  effects 

>  Motion  mechanical  systems  intrusive 

>  Flight  tests  costly,  time-consuming 

>  Computational  modeling  an  important 

element  for  advancing  fundamental 
understanding  and  engineering  prediction 


Unmanned  Combat  Air 
Vehicles  (UCAV) 


Simulation  provided  by  Mr  Ken  Wurtzler,  Cobalt  Solutions  LLC 


Massive  Separations/Vortic 

♦  Challenges  and  issues 


•  flow  fields  are  inherently  unsteady,  chaotic,  and  three- 
dimensional 

»  accuracy  is  crucial  at  high  angle  of  attack:  lift,  drag,  and 
moments 

»  complex  nature  of  massive  separation/vortical  flowfields 

•  defeats  conventional  turbulence  models 

•  higher  fidelity  computational  techniques  required 

•  flow  fields  are  described  by  the  Navier-Stokes 
equations 

»  analytical  solution  for  aircraft  not  possible 


ChoiflIETBffromputationalBodel 


♦  Direct  Numerical  Simulation  (DNS) 

•  solution  of  the  Navier-Stokes  equations  without  use  of  an  explicit 
turbulence  -  limited  to  low  Reynolds  numbers 

•  powerful  research  tool 

•  ready  for  full  aircraft  in  ~ 2080 


♦  Large  Eddy  Simulation  (LES) 

•  direct  resolution  of  the  large,  energy-containing  scales  of  the 
turbulent  flow,  model  only  the  small  eddies 

•  high  computational  cost  in  boundary  layers 

•  ready  for  full  aircraft  in  ~ 2045 

♦  Reynolds-average  Navier-Stokes  (RANS) 

•  model  the  entire  spectrum  of  turbulent  motions 

•  Highly  unreliable  performance  in  separated  flows 

•  ready  for  full  aircraft  today 


DES  method  combines  RANS  and  LES 


increase  in  empiricism 


Detached-Eddy  Simulation 

(DES) 

♦  Turbulence  modeling  approach  proposed  by  Spalart  et  al. 

(1997) 

•  Combines  Large  Eddy  Simulation,  and  Reynolds- 
Averaged  approaches 

•  Designed  to  provide  accurate  solutions  for  massively 
separated  flows 

•  Can  resolve  unsteady  flow  features 

»  Aero-acoustics,  aero -elasticity 

•  RANS  model  responsible  for  predicting  BL  growth 
and  separation  (NUMERICALLY  FEASIBLE) 

•  LES  model  responsible  for  prediction  of  unsteady 
flow  in  separated  region  (ACCURATE) 
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Simulation 
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❖  Scott  Morton  (PI),  Jim  Forsythe,  Tony  Mitchell 

♦>  AFOSR  project:  Aeroelasticity  predictions 
(PM:  Tom  Beutner,  John  Schmisseur) 

❖  AIAA  02-0587 
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Coarse  Grid 


X=500  mm 


X=800mm 


Real  Fine  Grid 


X=500  mm  X=600mm 


X=700mm  X=800mm 


Adaptive  Mesh  Refinement  Grid 
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Coarse  Grid  (1.2M  Cells) 

Medium  Grid  (2.7M  Cells) 

Fine  Grid  (6.7M  Cells) 
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Normal  Force  Power 


Spectral  Density  Analysis 
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♦  Data  of  Polhamus 

♦  Re=800,000 

♦  a=10° 

♦  Computations  made  on 
structured  and 
unstructured  grids  of 
various  domain  sizes  and 
grid  spacing 


Kyle  Squires  (PI),  Jim  Forsythe,  Philippe  Spalart 

AFOSR  project:  Spin  prediction 
(PM:  Tom  Beutner) 

DNS/LES  TV,  ERCOFTAC  Yol  8 


rounded-square  cross  section 
corner  radius  is  1/4  of  the 


6: 1  rectangular  ogive 
main  section  3.5b 


diameter 


endcap  0.5b 


DES 

pressure 


F-15E  at  65°  alpha 

rid  consists  of  5.9M  cells  (half  aircraft) 

•  Prisms  in  the  boundary  layer  (using  blacksmith) 

»  Conversions  to  prisms  saved  2M  cells 

•  Tetrahedrons  elsewhere 

•  Average  first  y+=0.7 

•  One  man-week  to  create 

•  Re=13.6xl06 

♦  2  days  to  compute  10,000  iterations  on  256 
processors  (tempest  -  MHPCC) 

♦  Time  step  and  grid  sensitivity  examined 

Jim  Forsythe  (PI),  Kyle  Squires,  Ken  Wurtzler,  Philippe  Spalart 

AFOSR  project:  Spin  prediction 
(PM:  Tom  Beutner) 
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Comparison  at  Station  4  (Spin  Coeff  =  0.2) 
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F-15E  at  65°langle 

♦grid  (full  aircraft):  6.46  x  106  cells  (generated  using 
VGRIDns) 

•  prisms  in  the  boundary  layers,  tetrahedra  elsewhere 

»  conversion  to  prisms  using  blacksmith 

•  average  first  y+  =  0.8 

•  Between  resolution  of  coarse  and  baseline  grids 

♦  timestep  =  0.02  (dimensionless  using  chord  length  and 
freestream  speed) 

+Re  =  13.6  x  106 ,  Mach  number  =  0.3 

♦  rotary  motion  about  centroid  along  freestream  velocity 
vector 


►  Bump  added  to 
nose  to  reproduce 
strong  yawing 
moment  seen  in 
flightiest 
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Detached  Eddy 
Simulation 


♦  importance  of  including  “LES  content”  in  the 
boundary  layer  prior  to  separation 


mbedded  LES  f of 
channel  flow 


Reseaf 
t 


•  flows  with  shallow  separation 


•  need  grid  densities  sufficient  to  sustain  eddy  content  near  the 

wall 


interface,  yj5  between  RANS  and  LES  regions  controlled  by  the  grid 


“super  buffer”  between  RANS  and  LES  velocity  profiles 
under-prediction  of  the  skin  friction  (Nikitin  et  al.  2000) 


Flow  Structure  near  RANS-LES  interface 


u  velocity  fluctuation,  DES,  129x129x65,  y+=250 
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Backscatter 


♦  stochastically  force  the  Navier-Stokes  equations  (Leith 
1990,  Mason  and  Thompson  1993,  Carati  et  al.  1994...) 
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f.  =  stochastic  force  distributed  about  RANS-LES  interface 


purely  random  or  scaling  using  the  eddy  viscosity,  strain  rate, 
and  timestep 


envelope  over  which  force  distributed 


Flow  Structure  near  RANS-LES  interface 


,  DES  with  backscatter,  129x129x65,  y+=250 
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Eddy  viscosity,  DES  with  backscatter,  129x129x65,  y+=250 
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Stress  comparison  for  DES  with  and  without  backscatter,  65x75x1 


stochastic  forcing  raises  resolved  stress,  lowers  modeled  stress 
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Future  Areas  of  Research 

Necessary 

♦  Embedded  LES  to  improve  simulation  of 
instabilities  generated  inside  the  boundary  layer 

•  Need  to  continue  the  research  outlined  above 

•  Apply  the  method  to  more  test  cases 

♦  Unsteady  experiments  of- 

•  Static  high  alpha  UCAV  configurations 

•  Pitch  and  roll  maneuver  tests  with  unsteady  data 
gathered 

•  Possibly  adopt  the  Boeing  1301  or  1303  as  a  standard 
configuration  for  several  groups  to  test 

•  High  accuracy  methods  applied  such  as  PIV,  LDV,  etc. 


♦  DES  has  been  examined  on  a  wide  range  of 
massively  separated  flows 

•  Moderate  to  greatly  increased  accuracy  over 
traditional  methods 

•  Capability  to  predict  unsteady  flows  at  flight  Re 

»  Crucial  for  high  alpha  maneuvering 
»  Crucial  for  aero-elasticity,  aero-acoustics 

•  Enough  confidence  built  to  encourage  engineering 
use 

♦  Several  areas  of  research  needed  to  apply  to 
super-maneuvering  UCAVs  with  confidence 


